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Recent studies have predicted extraordinary properties for transverse domain walls in cylindrical 
nanowires: zero depinning current, the absence of the Walker breakdown, and applications as domain 
wall oscillators. In order to reliably control the domain wall motion, it is important to understand 
how they interact with energy barriers, which may be engineered for example through modulations 
in the nanowire geometry (such as notches or extrusions) or as inhomogeneities in the material's 
crystal anisotropy. 

In this paper, we study the motion and depinning of transverse domain walls through potential 
barriers in ferromagnetic cylindrical nanowires. We use (i) magnetic fields and (ii) spin-polarized 
currents to drive the domain walls along the wire. Barriers are modelled as a section of the nanowire 
which exhibits a uniaxial crystal anisotropy where the anisotropy easy axis and the wire axis enclose 
a variable angle 9p. Using (i) magnetic fields, we find that the minimum and the maximum fields 
required to push the domain wall through the barrier differ by 30%. On the contrary, using (ii) spin- 
polarized currents, we find variations of a factor 130 between the minimum value of the depinning 
current density (observed for 8p = 0°, i.e. anisotropy axis pointing parallel to the wire axis) and 
the maximum value (for 8p — 90°, i.e. anisotropy axis perpendicular to the wire axis). 

We study the depinning current density as a function of the height of the energy barrier using 
numerical and analytical methods. We find that for an industry standard energy barrier of 40 kp,T, 
a depinning current density of about 5 fiA (corresponding to a current density of 6 x 10 10 A/m 2 in 
a nanowire of 10 nm diameter) is sufficient to depin the domain wall. 

We reveal and explain the mechanism that leads to these unusually low depinning currents. One 
requirement for this new depinning mechanism is for the domain wall to be able to rotate around its 
own axis. With the right barrier design, the spin torque transfer term is acting exactly against the 
damping in the micromagnetic system, and thus the low current density is sufficient to accumulate 
enough energy quickly. These key insights may be crucial in furthering the development of novel 
memory technologies, such as the racetrack memory, that can be controlled through low current 
densities. 



I. INTRODUCTION 

The current induced motion of domain walls (DWs) in 
ferromagnetic nanowires has been the subject of intense 
study in recent years. Being able to move and control 
DWs accurately is an important step towards the realisa- 
tion of devices such as the racetrack memory 1 . Two main 
challenges are the high critical current density, j'd, re- 
quired to depin the DW and initiate the DW motion 2 and 
the so-called Walker breakdown, a phenomenon which 
limits the DW velocity and is caused by deformations 
in the magnetization structure 3 . It has been recently 
observed 4 ' 5 that transverse domain walls (TDWs) in 
cylindrical nanowires are not subject to such limitations: 
j'd is zero, and the Walker breakdown does not occur. 
These DWs are able to propagate uniformly, without any 
deformation of their internal structure. Moreover, recent 
theoretical studies 6-8 suggest that TDWs may play an 
important role in the emerging research area of domain 
wall oscillators 9-13 (DWOs), providing an effective way 
to sustain magnetization dynamics with direct currents. 

The main obstacle for the experimental observation 
of these effects is that TDWs only occur in cylindri- 
cal nanowires with small diameter (below ~ 50 nm for 



Permalloy). In wires with greater diameter the demag- 
netizing field is strong enough to force the magnetiza- 
tion to follow the cylindrical surface of the wire, thus 
changing the internal structure of the DW, giving rise to 
a vortex DW 4 ' 14 . As a consequence, the experimental 
validation of the predicted properties of TDW requires 
the fabrication of cylindrical nanowires with diameter 
below 50 nm, which is challenging from an experimen- 
tal point of view. Recently, however, novel techniques 
have been developed for fabricating cylindrical nanowires 
with small diameter 15-17 . Ruitao et. al. 16 describe a pro- 
cedure to obtain single crystalline Permalloy nanowires 
with an average diameter of <~ 30 nm, through in situ fill- 
ing of the inner cavities of carbon nanotubes. Pitzschel 
et. al. 17 discuss the fabrication of arrays of Ni cylin- 
drical nanowires with diameter between 80 and 160 nm, 
using elcctrodeposition. Such arrays are technologically 
relevant, as they may represent a practical and conve- 
nient way to densely pack nanowires for memory storage 
applications. In particular, the authors discuss how the 
diameter of the nanowires can be modulated in order to 
introduce geometrical "artifacts" which act as barriers 
for the propagation of DWs. These barriers may be nec- 
essary in order to control the DW motion and make it 
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FIG. 1. (Color online) Sketch of the system. The arrows 
on the cylinder axis represent the magnetization. The trans- 
verse DW on the left is pushed through the barrier on the 
right. The barrier is modeled as a region where an additional 
uniaxial anisotropy pins the magnetization along the direction 
determined by Op. 



more predictable and reliable. It is then important for 
technological applications to study and understand the 
motion of TDWs through obstacles (which may be real- 
ized as inhomogeneities in the material or in the geometry 
in the nanowire). 

In this paper we present a micromagnetic study of 
the field-driven and current-driven motion of transverse 
domain walls through potential barriers in cylindrical 
nanowires. We consider two kinds of barriers: cylindri- 
cally symmetric barriers, whose associated energy term 
is invariant for rotations of the magnetization around 
the nanowire axis, and cylindrically asymmetric barriers, 
whose energy term does depend on the rotational state 
of the magnetization. When passing through a cylindri- 
cally symmetric barrier the TDW is "rotationally-free" , 
i.e. it can rotate freely around the nanowire axis. On 
the contrary, when passing through a cylindrically asym- 
metric barrier the TDW is "rotationally-bound" , as the 
barrier introduces an energy penalty for rotations of the 
DW. We find that the current density required to push 
the DW through the barrier is radically different in the 
two cases. In particular, we report that current densities 
of 10 10 - 10 11 A/ni can be used to depin rotationally- 
free DWs from barriers which are thermally "stable" 
at room temperature (the associated energy is greater 
than 40 k^T at T = 300 K). In contrast, depinning a 
rotationally-bound DW requires current densities which 
are higher by a factor ~ 130. 



II. THE SYSTEM 

In order to study the field-driven and current-driven 
motion of a DW through a potential barrier, we employ 
the system schematically shown in Fig. 1, a long cylin- 
drical ferromagnetic nanowire. The potential barrier is 
modeled as a thin region of the wire of length Lq where a 
uniaxial anisotropy pins the magnetization along a given 
direction a with a given strength K\ . In this first part of 
the paper, we fix the value of the anisotropy constant K\ 
and study how the motion of the DW depends on 6p, the 
angle between the unit vector a and the nanowire axis. 



The nanowire is subdivided into three regions: one 
longer region on the left with length Ll and one shorter 
region on the right with length L R , separated by the 
barrier, a thin layer of width Lq. In the simulation we 
"place" a DW on the left of the barrier and try to push 
it to the right using increasing applied fields (Sec. Ill) or 
current densities (Sec. IV). 

The dynamics of the magnetization, M, is computed 
using the Landau-Lifshitz equation, extended with two 
additional terms in order to model spin transfer torque 
effects 18 : 



d t M = — 7 M x H + — M x d t M 

M s 

+ vd x M-^rMxd x M (1) 



In this equation, M s = \\M\\ is the saturation magnetiza- 
tion, H is the effective magnetic field, 7 is the gyromag- 
netic ratio, a is the damping parameter. The current 
density, ja, pp , is applied along the x direction and enters 



-, where 



the model through the parameter v = eM ,( 1+ g2) : 
P is the degree of polarization of the spin current, [Ib 
is the Bohr magneton, e the absolute value of the elec- 
tron charge, £ the non-adiabatic parameter. Note that in 
this paper we reason in terms of fully polarized current 
densities, j = Pj app , rather than in terms of the actual 
applied current densities, j app . 

The effective field is calculated as H = i? C xch + 
-ffdcmag + H app + ffpin, where H cxch = V 2 m is the 

exchange field, -ffdcmag is the magnetostatic field, H app 
is the applied field and -ff p i n = (2K±m ■ a/fioM s ) a is the 
magnetic anisotropy corresponding to the pinning poten- 
tial U p i n = —K\ (m- a) 2 . K\ is the only parameter which 
varies in space: it is zero outside the barrier region, and 
is set to K\ = 0.5 x 10 6 J/m 3 inside it. The other param- 
eters are homogeneous throughout the whole nanowire 
and are set to M s = 0.86x 10 6 A/m, A = 1.3x 10~ n A/m, 
7 = 2.21 x 10 5 m/(As), a = 0.01 and £ = 0.01, which 
are typical values for Permalloy. Concerning the geom- 
etry of the cylindrical nanowire we take Ll = 148.5 nm, 
Xr = 48.5 nm, Lq = 3nm, so that the total length is 
L = Ll + Lq + Ljx = 200 nm. The diameter is chosen 
as D = 10 nm. Contributions from the Oersted field and 
the Joule heating are negligible for nanowires of small 
radius 7 and are thus not included in the model. 

We note that, when 6 P > 0°, the pinning field H p [ n 
tries to pull the magnetization out of the axis of the 
nanowire. The magnetostatic field, however, opposes this 
and — for the value of Ki considered in this paper - 
manages to keep the magnetization along the axis of the 
nanowire, even inside the potential barrier. 

Eq. (1) is discretized over a finite element mesh and 
the simulations are carried out using a version of the 
micromagnetic simulation package Nmag 19,20 , extended 
to model spatially varying magnetic anisotropies. 
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III. FIELD-DRIVEN DW MOTION 

We carry out a number of simulations to determine 
the critical field, i? C rit, which must be applied (in the 
direction of the nanowire axis) in order to push the DW 
through the potential barrier. In particular, we fix the 
pinning strength to K\ = 0.5 x I0 6 J/m 3 and perform one 
simulation for each value of the pinning angle # P , which 
is changed from zero (pinning along the nanowire axis) 
to 90° (pinning orthogonal to the nanowire axis) in steps 
of 5°. 

Each simulation is carried out in two parts. In the first 
part, a preliminary relaxation computation is carried out 
and the system settles into a metastable state where a 
DW is located on the left of the barrier. In the second 
part, the applied field is increased gradually until the 
DW passes through the barrier. For the preliminary re- 
laxation we set the initial magnetization in the following 
way, 



M(x,y,z) = M s (sin0(x), cos 0(b) , 0), 



(2) 



where 9(x) = | min (1.0, max (— 1.0, (x — xq)/w)) and 
xq = — 30 nm, w = 20 nm (the reference system is chosen 
as in Fig. 1). For this choice of parameters, the magne- 
tization relaxes into a tail-to-tail DW, located on the left 
of the nanowire, similarly to what is schematically shown 
in Fig. 1. The components of the magnetization along 
the axis of the nanowire before and after the relaxation 
computation are shown in Fig. 3. The relaxation is car- 
ried out using an artificially high value for the damping, 
a = 0.5, to speed up the computation, and is done in the 
presence of a small magnetic field, Ho = 10 3 A/m, which 
is applied along the negative x direction and is used to 
push the DW towards the barrier. We want indeed to 
ensure that the DW settles into an equilibrium position 
against the barrier (e.g the center of the barrier) before 
starting to "depin" it. The convergence criterion used to 
stop the simulation is the following: 



max ■ 



dM, 



df 



< e 



(3) 



where the index i runs over all the sites of the finite 
element mesh and e = 1 deg/ns. 

The procedure described above brings the system to 
a state which is different for different values of the pin- 
ning angle, 0p. In particular, when P is zero the DW 
approaches the barrier from the left without penetrating 
it (Fig. 2(a).), as this would require the magnetization 
in the region to move out of the pinning direction. On 
the other hand, when P is 90° the DW "falls" into the 
center of the barrier allowing the magnetization to align 
along the pinning direction, with a consequent reduction 
in energy (see Fig. 2(b)). Depending on 9p we then have 
two different scenarios, where the barrier either repels 
the DW (as for P = 0°) or attracts it (as for P =90°). 

The magnetization configuration obtained from the 
preliminary relaxation computation is used as the initial 
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FIG. 2. (Color online) The relaxed magnetization configura- 
tions for Op = 0° (a) and 0v = 90° (b). The sketches on the 
corners of the two figures show that the DW moves into the 
barrier to allow the magnetization to align along the pinning 
direction in the barrier. 
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FIG. 3. (Color online) The components of the normalized 
magnetization, m = M/M s> along the axis of the nanowire 
before (dashed lines) and after (solid lines) the preliminary 
relaxation computation. The blue curves (light in greyscale) 
show the x (axial) component of m, while the dark red curves 
(dark in greyscale) show the component of m orthogonal to 
the axis. 



state for the second part of the simulation, in which the 
damping is set to the realistic value, a = 0.01, and the 
applied field is gradually increased until the DW passes 
completely through the potential barrier. In particular, 
we increase the field in steps of AH = 10 2 A/m and for 
each value of the applied field we relax the system, i.e. 
perform a time integration of Eq. (1) until the conver- 
gence criterion (3) is met. To determine when the DW 
passes through the barrier and stop the simulation, we 
analyze the values of the spatially averaged normalized 
magnetization (m), where in = M/M s . The x compo- 
nent of this vector, {m x ), provides a good indication of 
the position of the DW, while the y and z components 
provide information about its state of rotation. In partic- 
ular, if we assume the DW to be mirror-symmetric along 
the x axis, we get (m x ) ~ (Lp — £l — 2xjjw)/L, where 
Xr,w is the position of the DW, while L is the total length 
of the nanowire, L = Ll + Lq + Lp. In the simulations 
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FIG. 4. (Color online) H CI i t , the field required in order to 
push the DW throughout the potential barrier, is plotted as 
a function of the pinning angle, Op. 
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FIG. 5. (Color online) jait, the current required in order to 
push the DW through the potential barrier, is plotted as a 
function of the pinning angle, Op. 



we assume the DW to pass through the potential barrier 
when (m x ) < —0.75, which corresponds, following the 
formula above, to a position of idw = 25 nm. 

The results of the simulations are collected in Fig. 4. 
The field necessary to move the DW through the barrier 
varies between (68.05 ± 0.05) and (106.55 ± 0.05) kA/m. 



IV. CURRENT-DRIVEN DW MOTION 

In order to determine the current necessary to push the 
DW through the barrier we follow a procedure analogous 
to the one described in the previous section, with two 
main differences. First, to push the DW we use a current 
rather than a field. Second, we abandon the convergence 
criterion in Eq. (3). This convergence criterion works 
well in the field-driven case because the torque is guar- 
anteed to decrease in time due to the damping term. In 
the current-driven case, however, the spin transfer torque 
can oppose the damping and bring the system to a steady 



state where the magnetization dynamics is sustained in 
time. In particular, it may happen that the DW reaches 
a state where it does not advance anymore along the wire 
but keeps rotating around its axis 7 . In this situation, the 
torque never falls below the given threshold, e, and Eq. 
(3) is never satisfied. We then use a convergence criterion 
which stops the simulation when the DW translational 
velocity (along x) falls below a certain threshold. For 
convenience, we formulate it in terms of {m x ), as follows: 



(m x (t + At)) - (m x (t)) 



At 



< e, 



(4) 



with e = 10~ 5 ns _1 . The convergence criterion is checked 
every At = 100 ps and the simulation is stopped when it 
is satisfied for five consecutive times. We note that Eq. 
(4) may be expressed in terms of the DW velocity using 
the relation (m x ) — (L-r — Ll — 2xr)w)/L, 



"DW 



xnw(t + At) - a; DW (Ai) 



At 



L 



(•5) 



which reduces, in our specific case, to i>dw < 
10~ 3 nm/ns. We also pin the magnetisation at the left 
and right ends of the simulated wire to reduce finite-size 
effects (such as the nucleation of DW, which may occur 
at current densities above ~ 10 12 A/m 2 ). 

As in the field-driven case, we carry out one simula- 
tion for each chosen value of 9p and each simulation is 
subdivided into a preliminary part, where the DW is pre- 
pared in its initial state on the left of the barrier, and a 
main part, where the current density is increased until 
the DW passes through the barrier. In the preliminary 
simulation we use a current density jo = 10 10 A/m 2 to 
push the DW from the left to the right against the poten- 
tial barrier. The current is directed along the negative 
x direction. We use a high damping, a = 0.5, simi- 
larly to the field-driven case. In the main part of the 
simulation, we restore the damping to the more realistic 
value a = 0.01 and increase the current density in steps 
of Aj — 10 10 A/m 2 . The results of the simulations are 
shown in Fig. 5. 

We see that the current density required to push the 
DW through the potential barrier, j cr ; t , increases rapidly 
as the pinning direction moves out of the axis of the 
cylinder: for P = 0° we get j clit = (0.0239 ± 0.0001) x 
10 12 A/m 2 (this value was determined with a separate 
simulation to get improved accuracy) , while for 9 P = 90° 
we get j cri t = (3.185 ±0.005) x 10 12 A/m 2 . We conclude 
that the critical current j cr ; t varies roughly by a factor 
130 depending on the direction along which the barrier 
pins the magnetization. This is a quite remarkable fea- 
ture considering what happens in the field-driven case, 
where H cl i t does actually decrease for pinning angles ap- 
proaching 90°. 

To understand why the value of Op has such a strong in- 
fluence on the critical current j cr n we have to analyze the 
dynamics of the magnetization. Fig. 6 shows the time 
evolution of the components of the spatially averaged 
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FIG. 6. (Color online) Comparison of the current driven 
dynamics of the spatially averaged magnetization, (fn), for 
9 P = 0° (top) and P = 20° (bottom). The inset in the bot- 
tom plot shows the oscillatory motion of the DW induced by 
the pinning in the barrier. The discontinuities in (m,) corre- 
spond to changes in the applied current density (see Fig. 7) 
following the procedure described in the text. Two different 
scales are used for the time axis of the bottom plot in order to 
give clear information on the dynamics during the two differ- 
ent phases of the switching process. Note that the dynamics 
of the preliminary relaxations is omitted in these figures. 
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FIG. 7. (Color online) The applied current density for the 
two simulations of Fig. 6 and Fig. 8. The vertical lines in the 
plot mark a change in the scale of the time axis. 



magnetization, (fn), during the simulations for Op = 0° 
(top) and 6 P = 20° (bottom). When 9 P = 0°, the DW 
rotates around the nanowire axis, as shown by the y and 
z components of (fn). At the same time, the x com- 
ponent shows that the DW gradually compresses against 
the barrier as the current density, j, is increased from the 
value 10 10 A/m 2 , to 2 x 10 10 A/m 2 and 3 x 10 10 A/m 2 at 
times ns, 54.9 ns and 89.4 ns, respectively (dashed line 
in Fig. 7) . The effect of the spin transfer torque is con- 
stantly accumulating in time and is stored in form of 
barrier penetration and DW compression. 

In the case Op = 20°, the potential barrier pins the 
magnetization out of the direction of the nanowire axis, 
thus breaking the rotational symmetry of the system (i.e. 
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FIG. 8. (Color online) Evolution of the energy for the two 
cases 8p = 0° (dashed line) and 8p = 20° (solid line) con- 
sidered in Fig. 6 and 7. Units of k B T = 4.14 x 1CF 21 J 
(T = 300 K) are used for the energy. The time axis uses 
different scales as in Fig. 7. 

the energy of the system is not invariant for rotations of 
the magnetization). As a consequence, the DW cannot 
freely rotate around the axis of the nanowire, but rather 
precesses narrowly around its equilibrium configuration. 
The inset at the bottom of Fig. 6 confirms that (m x ) 
oscillates weakly and does not vary significantly during 
the application of increasing current densities (solid line 
in Fig. 7). Only at time 791.0 ns, when the current 
density reaches the critical value, j cr ;t = 24 x 10 10 A/m 2 , 
the DW starts to penetrate significantly into the barrier. 
In then takes about 2 ns for the DW to pass completely 
through it. 

Fig. 8 shows the micromagnetic energy, AE, in the two 
cases Op = 0° (dashed line) and 9p = 20° (solid line). 
AE(t) is computed as E(t) — E(0), where E(t) is the 
energy of the system at time t. The height of the energy 
barrier, AE max , is almost the same in the two cases. In 
particular, for P = 0°, A£ max = 39.24 k B T while for 
Op = 20°, Ai? max = 39.16 AsbT. It is interesting to notice 
that in the case of a rotationally symmetric barrier, the 
energy increases gradually as a response to the increase 
of current density. On the contrary in the case of an 
asymmetric barrier, the energy changes very weakly with 
j. In particular, the energy changes less than 1.2 k^T 
when j is increased from 10 10 A/m 2 to 23 x 10 10 A/m 2 . 
At time t = 791 ns, however, the current is increased to 
the value j cr it = 24 x 10 10 A/m 2 and it takes only 5 ns 
for the DW to pass through the remaining part of the 
barrier. 



V. DISCUSSION 

To understand why the rotationally symmetric barrier 
can be penetrated much more easily than the asymmetric 
one, we start from Eq. (1). We express the equation in 
the form where the time derivative of fh appears only 
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on the left hand side (this requires a substitution of the 
equation in itself) and both of its sides have been divided 
by lM s : 

— d t rh — —rh x H — am x (rh x H) 

i 



+ (l + aO-d x m 
7 

— (£ — a) — rh x d x rh, 



(6) 



where 7' = 7/ ( 1 + a 2 ) and m = M /M s is a unit 
vector. The right hand side of this equation receives 
the contributions of four terms. They are, respectively, 
the precession term, the damping term, the adiabatic 
Spin Transfer Torque (STT) and the non-adiabatic STT. 
The strongest of these terms is typically the preces- 
sion term, which receives contributions of the order of 
10 6 A/m (from the applied field, the demagnetizing field, 
the magneto-crystalline anisotropy) or higher (from the 
exchange field). The damping term has lower magnitude, 
due to the prefactor a. For Permalloy, a ~ 0.01 and the 
damping term is then 100 times weaker than the preces- 
sion term. Despite this, the damping term has a very 
important role for the physics of the system and has sub- 
stantial impact on its dynamics. The damping is very 
effective because it keeps changing the magnetization to- 
wards the direction that yields the highest reduction in 
energy. 21 Its effect, hence, accumulates in time and brings 
the system to a local energy minimum. 

On the other hand, the first STT term (the second STT 
term is smaller of a factor £ — a and will be neglected in 
this analysis) does not have this property: while it can 
have magnitude similar to the damping term (for Permal- 
loy and j = 10 12 A/m 2 , the adiabatic term can bring 
contributions of the order of 10 4 A/m in Eq. (6)), it does 
not always point in the direction which is best suited to 
pump energy into the system. Often, this term acts in the 
direction opposite to the much stronger precession term, 
with the result that the STT pumps energy in and out 
in an alternate fashion, i.e. there is no "accumulation" 
of the effect. 

Fig. 9(a) provides an example of this phenomenon. 
The plot shows the initial part of the trajectory of 
the spatially averaged normalized magnetization (rh) = 
(M/M s ) in the plane yz for the case 9 P = 20° of Fig. 6. 
The points of the trajectory where the energy of the sys- 
tem is increasing, dfE > 0, are coloured in red; the points 
where it is decreasing, d t E < 0, are coloured in green. 
The line thickness is proportional to |c?t-E|: a thick red 
line corresponds to maximum energy increase, while a 
thick green line to maximum energy decrease. The plot 
shows that once the current is applied, the energy starts 
increasing while the magnetization moves in the nega- 
tive y direction. At this point, however, the pinning is 
strong enough to oppose the current-driven dynamics and 
pulls the magnetization back, giving rise to a motion spi- 
ralling towards a new equilibrium configuration. During 
the precession, the STT pumps energy in and out, in an 
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FIG. 9. (Color online) Current driven trajectories of the aver- 
age magnetization, (rh), projected on the yz-plane (orthogo- 
nal to the nanowire axis). The dynamics is in response to the 
application of a current density j — 2 x 10 10 A/m 2 in the case 
of asymmetric barrier, 6p = 20° (a) and symmetric barrier, 
8p = 0° (b). The line thickness is proportional to the mag- 
nitude of the energy which is being pumped in (red in color, 
dark in greyscale) or out (green in color, light in greyscale) at 
the considered point. The two black dots indicate where the 
trajectories begin. Note that side (a) shows a few precessions 
of {m x ), while side (b) shows only one precession. This is just 
to avoid the overlapping of the line. 



alternate fashion. Overall, the current produces changes 
in rh of the order of 4 x 10~ 4 radians. 

A completely different situation is shown in Fig. 9(b), 
which shows again the dynamics of (rh), but for a DW 
moving through the rotationally symmetric barrier (case 
dp = 0° in Fig. 6). The plot shows that the energy is al- 
ways pumped in (and never pumped out). Moreover, the 
rate with which this happens (line thickness) increases in 
time. This is due to the compression of the DW, which 
causes a stronger STT effect (higher values of d x rh in Eq. 
(6)). 

The analysis we have carried out demonstrates that it 
is of fundamental importance to find systems and situ- 
ations which allow the STT to coherently pump energy 
into the system. The effectiveness of the STT in the 
case of rotationally symmetric and asymmetric barriers 
depends on whether or not the STT points along the op- 
timal direction for energy pumping (i.e. opposite to the 
damping) . 



VI. ENERGY BARRIER 

It may be useful at this point to give a more quanti- 
tative indication on the height of the barriers which can 
be penetrated thanks to the mechanism described in this 
paper. In particular, we would like to find out, given a 
rotationally symmetric barrier with a given height AE, 
what the required critical current density is in order for 
the DW to pass through it. 
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FIG. 10. (Color online) Study of the critical current, j cr it 
required to push the DW through a rotationally symmetric 
barrier, as a function of the barrier height (AE). The top 
x axis shows the values of the anisotropy constant K\ corre- 
sponding to the values of AE shown on the bottom axis. 



We perform simulations similar to the ones described 
previously, but here we fix 6p = 0° and change the value 
for the pinning anisotropy Ki from 0.1 x 10 6 J/m 3 to 
1.1 x 10 6 J/m 3 in steps of 0.1 x 10 6 J/m 3 . For every sim- 
ulation we also calculate the energy height AE as the 
difference AE — E max — E(0), where E(0) is the en- 
ergy of the system at the beginning of the simulation, 
while E max is the maximum value reached throughout 
the simulation. The values of the critical current, j C rit, 
are plotted in Fig. 10 as function of the barrier height, 
AE. The corresponding values for K\ are also shown at 
the top of the plot. We note that there is a linear rela- 
tionship between all the three quantities involved: jcrit, 
AE and K\ . In particular, AE and K\ (bottom and top 
x axis) are linearly proportional and follow the relation 
AE(Ki) = V K-y where V is a parameter with dimen- 
sions of a volume and V — (326.4 ± 0.4) nm 3 (from a 
least squares fit). This means that — as one may expect 
- a doubled value for K\ corresponds to twice-as-high 
energy barrier, AE. 

The linear relation between j cr ;t and AE is somewhat 
more difficult to understand intuitively. Fig. 10 shows, 
however, that the data can be well fitted (see dashed 
curve) against the function, 



j clit (AE) = BAE. 



(7) 



In particular, we get B = (1.498± 0.006) x 10 29 A/(m 2 J). 
Or — in terms of thermal energy — it is necessary to 
apply a fully polarized current density j = Pj apP of 2.48 x 
10 10 A/m 2 in order to push the DW through a potential 
barrier with height 40 k^T. Assuming a spin polarization 
P = 0.4, we find j app = 6.20 x 10 10 A/m 2 . 

We now multiply both sides of equation (7) by the cross 
sectional area of the nanowire, irR 2 : 



I crit (AE)=BAE, 



(8) 



which relates the critical current / C rit — irR jcmt 

to AE. 

We have B = ttR 2 B = (1.177 ± 0.005) x 10 13 A/J. 



VII. ANALYTICAL CONSIDERATIONS 

In this section we derive a simple analytical model with 
the aim of understanding the linear proportionality be- 
tween the current density, j cl n, required to push a DW 
through a potential barrier, and the barrier height, AE 
(as seen in Fig. 10). We start from the analytical model 
developed in Ref. 7 (the analytical derivation below is 
heavily based on this previous work). In particular, we 
consider the situation where the pinning potential is in- 
finitely high and the DW does not penetrate into the 
barrier, but rather compresses against it, thus accumu- 
lating a certain amount of exchange/compression energy. 
It is reasonable to expect that if the applied current can 
accumulate an amount of "compression energy" , AE, in 
the case of an infinite barrier, then it may also be able to 
push the DW through a finite energy barrier of the same 
height. We stress that in this context we are just trying 
to get to a qualitative understanding of the phenomenon, 
rather than a quantitative analytical model. 

First, we write down the energy of the DW, using 
the one dimensional model, the coordinate system and 
nomenclature of Ref. 7, 



U = A 



[(d x 8) 2 + sin 2 6 (d x <b) 2 ] dx. 



U has the units of energy per cross-sectional area. A is 
the exchange coupling constant, L the DW width, 9 and 4> 
are spherical coordinates for the magnetization. In this 
expression, we consider only the exchange energy and 
neglect the contribution of the magnetostatic field. The 
term d x <j> is usually smaller 7 than d x 9, and we neglect it: 

U=~ [ {d u 6) 2 du = -y f d u 0d9, 
L Jo L Jo 

where we also changed the variable of integration twice, 
from itou = x/L and from u to 6. In the high current 



regime 



d u 



3V 



sin 2 , and, 



U 



3AV 
La 



sin — d9 



6AV 
La 



Since, V = ^ w, we have: 



U = 



6AV 



La 



e ja 



(9) 



For the low current regime the DW deformation is negli- 
gible and we get the energy of a relaxed transverse wall of 
length L, Uq = Air 2 /L. In summary, the present analysis 
is valid only in the high current regime, i.e. for currents 



j = Pj apP > jo 



2e"f aA 



MoMs L 
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In this regime, the DW compresses accumulating an ex- 
change energy U per unit of cross sectional area, 

U=\u»^, forj»j . 

71" JO 

The energy pumped in by the current is AU = U — Uo 1 

AU = U ( 4 - - l), for J ». j . (10) 
V 71 " Jo J 

In our system we have jo = 0.026 x 10 12 A/m 2 (we have 
assumed a DW length L — 30 nm). We expect Eq. (10) 
to hold only when the current exceeds this value. When 
j ~ jo, on the other hand, we may still have important 
STT effects, but getting to an analytical expression would 
require to solve exactly the system of equations (14) in 
Ref. 7. We can now compare the value obtained for B 
(see Eq. (8)) from the fit, B = (1.177±0.005) x 10 13 A/J, 
with B = 1.013 x 10 13 A/J, as obtained from the equation 
below (derived from Eq. (9)), 

B = Cia 

Note that B does not depend on the considered geometry, 
meaning that the total depinning current, J cr it, does not 
depend on the size of the nanowire. In larger nanowires 
then lower current densities can be used to overcome bar- 
riers of the same height in energy. For example, the cur- 
rent density, j app = 6.20 x 10 10 A/m 2 , required to push 
a transverse domain wall through a barrier of 40 k^T in 
a nanowire with diameter 10 nm, should be reduced by 
a factor 9 in nanowires with diameter 30 nm. The op- 
timization of the system geometry necessary to achieve 
lower current densities and lower switching times are left 
to future investigations. 



In summary, Eq. (10) shows that the amount of 
exchange energy the DW accumulates by compressing 
against the infinite potential barrier is proportional to 
the applied current density. Fig. 10 is then easily un- 
derstood, assuming the linear behaviour is maintained in 
the case of finite barrier. 



VIII. SUMMARY 

We have studied a transverse DW in a cylindrical 
nanowire and its motion through a barrier, modeled as a 
pinning potential on the magnetization. We determined 
the critical fields and current densities required to push 
the DW through the barrier for various directions of the 
pinning potential. We found that the critical applied field 
decreases as the pinning direction gets orthogonal to the 
nanowire axis. On the contrary, the critical current den- 
sity increases by more than a factor 130 when the pinning 
direction gets orthogonal to the nanowire axis, meaning 
that the DW can penetrate the barrier much more easily 
when the pinning potential is aligned along the axis of 
the wire, rather than being orthogonal to it. 

This study gives further insights into the extraordinary 
properties of transverse DWs in cylindrical nanowires and 
motivates experimental investigations on these systems. 
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